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Abstract. This is the first paper presenting our long-term project aimed at studying the nature of bulges through the analysis 
of their stellar population gradients. We present deep spectroscopic observations along the minor axis and the data reduction for 
a sample of 32 bulges of edge-on spiral galaxies. We explain in detail our procedures to measure their dynamical parameters 
(rotation curves and velocity dispersion profiles) and line-strength indices, including the conversion to the Lick/IDS system. 
Tables giving the values of the dynamical parameters and line-strength indices at each galactocentric radius are presented (in 
electronic form) for each galaxy of the sample. The derived line-strength gradients from this dataset will be analyzed in a 
forthcoming paper to set constraints on the different scenarios for the formation of the bulges. 
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1. Introduction 

Bulges of spiral galaxies are cornerstones for constraining the- 
ories of galaxy formation. Located at the centres of spiral 
galaxies, they hold the signature of the sequence of formation 
of the different sub-systems building a spiral galaxy: halo, disc, 
and bulge. While the mass of disks of spiral galaxies in the local 
universe is approximately constant among spirals of different 
Hubble types (e.g., Arimoto & Jablonka 1991), the prominence 
of bulges varies widely along the Hubble sequence. Hence, 
bulges are likely to hold important clues to our understanding 
of galaxy evolution and the nature of the Hubble sequence. 

Formation scenarios for bulges can be divided into 
two categories. One of them states that bulges are formed 
in a way similar to that of low- to intermediate-mass el- 
liptical galaxies, based on the strong similarities between 
global properties of many bulges and of elliptical galaxies. 
Bulges and (small) elliptic als populate the same location 
in the Fundamental Plane (Bender, Burst ein. & Faberl 1 1992 : 



Falcon-Barroso. Peletier. & Balcells 
a continuous sequence in the V max 



2002), and they form 
/ a vs. ellipticity di - 



Benderetal. 1992) 



agram (being supported by rotation; 
Furthermore, spectroscopic studies of the central part of bulges 
have shown that their mass-metallicity relation (when derived 



* Based on observations collected at the Roque de los Muchachos 
Observatory and at the European Southern Observatory, proposal 
numbers 58A-0192(A), 59A-0774(A) and 61A-0326(A) 



from a elements, e.g. the Mg 2 index) is consistent with 
that of elli ptical galaxies ( Jablonka. Martin. & Arimotol 



19961: 



1996; 



Jdiart. de Freitas Pacheco. & Costal 
Moo rthv & Holtzmanl 12006), 3ut see also Prugniel, Maubon 
& Simien (2001), and Proctor & Sansom (2002). We refer the 
reader to the introduction in Jablonka, Gorgas, & Goudfrooij 
(2007, hereafter Paper II) for a discussion on previous 
spectroscopic studies of bulges. 

The other popular type of formation scenario for bulges is 
the "secular evolution" scenario in which bulges form from 
disk material through redistribution of angular momentum. 
This scenario involves stellar bars driving gas from the disk into 
the c entral region of the galaxy, thus triggering star formation 
(e.g.JPfenniger & Normanlll99C : Norman. Sellwood. & Hasan 
19961: iKorm endv & KennicuttOIo04j). If enough mass is ac- 
creted, the bar itself will dissolve and the orbits of stars in- 
volved in the process will yield a bulge-like spatial distribu- 
tion. Galaxies would thus evolve from late to earlier types 
along the Hubble sequence. Observational support for secu- 
lar evolution has been provided mainly for late-type spiral 
galaxies (~ Sb and later), inclu ding correlations between the 
scale l engths of bulges and di sks (Courteau, de J ong. & Broeils 



19961 iMcArthur. Courteau. & Holtzmanl 120031) and the radial 
light distribution of bulges of late-type spirals being closer 
to exponential than t o classical de Vauc ouleurs' i? 1 / 4 profiles 
dBalcells et al.ll2003b Ide Jong et al.ll2004 . 
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Stellar population studies should be able to place impor- 
tant constraints on the formation mechanisms for bulges. If 
bulges form mainly through dissipative collapse without signif- 
icant subsequent merging, simulations indicate that the slope 
of radial metallicity gradients should steepen with increasing 
galaxy mass and luminosity (lCarlbergill984l) . Simulations of 
dissipative merging of gas-rich disks has also been shown to 
prod uce metallicity gradient s that steepen with increasing mass 
(e.g jBekki & Shiovall 1999b . Conversely, the impact of secular 
evolution on population gradients is expected to be rather dif- 
ferent. Simulations including gas and effects of star formation 
show that disk gas fueled to the central regions renders a cen- 
tral region that is younger and more metal rich than the outer 
regions, whereas metallicity gradients outside the central re- 
gion flatten out signi ficantly ~ 1 Gyr after formation of a bar, 
both for gas and stars ( Friedl iT~Benz. & KennicutJl994tlF^iedlil 
19981) . Observational evidence for this effect exists, albeit so 
far only for the gas component in spiral disks. Several studies 
have shown that global radial gradients of the gas metallicity 
in barred spirals are shallower than gradients in unbarred spi- 
rals of the same Hubble type ( Vila-Costas & Edmunds! 1 1992t 
Zaritskv. Kennicutt, & Huchralll994HMartin & Rovll 19941) . 

In terms of stellar population gradients in bulges, past 
studies mainly us e d bro ad-b and imaging. The studies of 
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Balcells & Peletierl d 1994b and IPeletier et ail d 1999b used op- 
tical and near-infrared colors and found that color gradients 
of luminous bulges increased with luminosity whereas faint 
bulges showed stronger gradients than expected from the trend 
observed for the luminous bulges. However, information in col- 
ors is generally rather degenerate in age, metallicity, and/or ex- 
tinction. 

Spectroscopic studies allow the measurement of line 
strength indices that are insensitive to dust extinction, and 
allow a cleaner separation of age and metallicity of a stel- 
lar population. While line strength gradients have been mea- 
sured for many elliptical and SO galaxies using the Li c k/IDS 
system of indices (e.g., ICarollo. Danziger. & BusorJ I l993[ 
Gonzalez & Gorgas 1996; Fisher, Franx, & Illingworth 1996; 
Meh lert et al.ll2003TlKuntschner et al.ll2006b . such data are rel- 
atively hard to acquire for bulges of spirals at high enough 
signal-to-noise ratios, given the often low surface brightness 
of bulges and the difficulty of disentangling the observed spec- 
tra into bulge and disk components. This is especially difficult 
for late-type spirals in face-on or moderately inclined configu- 
rations. 

With this in mind, we embarked on an extensive spectro- 
scopic survey of 32 spiral galaxies in an edge-on configuration 
to avoid contamination by disk light, and with Hubble types 
ranging from SO to Sc. Early results from subs ets of the data 



acquir ed during this su rvey w ere published in Goudf rooii et al 
dl999b . |Jablonka et al.ld2002l) . and lGorgas et al. Id2003b . In this 
first paper on the final and comprehensive results of this sur- 
vey, we present the extensive data reduction and analysis pro- 
cedures leading to the final line strength measurements. The 
galaxy sample is presented in Section [2] Section [3] describes 
the observations, and Section [4] provides details on the data 
reduction procedures. Section [5] contains a description of the 
determination of dynamical parameters (radial velocities and 
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Fig. 1. Histogram showing the distribution in Hubble types of 
the galaxy sample. 



velocity dispersions as a function of bulge radius), and the 
line strength index measurements themselves are described in 
Section [6] Finally, Section |7]contains a summary of this work. 

2. The sample 

We selected a sample of 32 genuine (or close to) edge-on 
spiral galaxies. Galaxies in the northern hemisphere were se- 
lected from the Uppsala General Catalogue (Nilson 1973), 
while southern galaxies were selected from the ESO/Uppsala 
catalogue (Lauberts 1982). As shown in Figure Q] the galax- 
ies span a large range in Hubble types (from SO to Sc). Given 
the edge-on configuration of the galaxies, precise morphologi- 
cal classification is difficult. Due to this we consider a Hubble 
T-type uncertainty of ± 1, 

which is also a fair representation of the catalogue-to- 
catalogue variations for a given galaxy. 

All the galaxies are nearby objects, with radial velocities 
between 550 and 6200 km s _1 , and were chosen to cover a 
range in bulge luminosity (—17.5 > My > —20.5). Table [T] 
lists the galaxy sample and some relevant parameters, including 
the corresponding bulge types (from Lutticke et al. 2000) and 
the local environment, which also varies significantly among 
the sample galaxies. A few galaxies have previous photomet- 
ric and spectroscopic observations from other authors (see 
Section^ and were included for comparison purposes. 

3. Spectroscopic observations 

Spectroscopic observations of the galaxy sample were obtained 
during five observing runs at three different telescopes: 

The northern galaxies were observed at the 2.5-m Isaac 
Newton Telescope (INT) of the Isaac Newton Group of tele- 
scopes on the island of La Palma (Spain). For the southern 
galaxies we used the 3.6-m ESO telescope and 3.6-m New 
Technology Telescope (NTT), both at ESO La Silla observa- 
tory (Chile). 
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Table 1. The galaxy sample. 





Type 


T 


J3t 


i 


Bulge Type a 


Environment 


run b 


PA C 


exp. d 


r cff 


Ar c ff 


V* 

v max 


NGC 522 


Sc 


4.1 


13.98 


90 


box 




2 


123 


8200 


3.5 


0.7 


196.1 


NGC 585 


Sa 


1.0 


14.06 


90 


close to box 




4.5 


86 


7200 


2.4 


0.7 


246.6 


NGC 678 


SBb 


3.0 


13.31 


90 


box 


cluster 


1 


168 


12000 


23.3 


1.4 


169.0 


NGC 891 


SAb 


3.0 


10.83 


88 


box 




1 


112 


10500 


13.2 


0.7 


212.0 


NGC 973 


Sbc 


3.2 


14.17 


90 


box 




2 


138 


13800 


3.4 


0.4 


270.0 


NGC 1032 


SO-a 


0.4 


12.74 


90 


elliptical 


companions? 


1 


158 


9000 


11.4 


1.1 


283.9 


NGC 1184 


SO-a 


0.0 


12.95 


90 


close to box 




2 


79 


8200 


11.3 


1.9 




NGC 1351 A 


SBbc 


4.5 


14.16 


80 


elliptical 


cluster 


5 


132 


5400 


5.8 


0.8 


90.7 


NGC 1886 


Sab 


3.9 


13.80 


90 


peanut 


isolated 


3 


60 


10800 


1.6 


0.3 


154.6 


NGC 3957 


SAO 


-1.0 


13.03 


90 


close to box 


cluster 


3 


173 


7200 


4.0 


0.8 




NGC 5084 


SO 


-1.8 


11.55 


90 


elliptical 


cluster 


3 


80 


7200 


38.1 


11.6 


310.3 


NGC 6010 


SO-a 


0.4 


12.95 


90 






1 


15 


13800 


4.2 


0.9 


148.8 


NGC 6829 


Sb 


3.0 


15.09 


90 






2 


121 


12800 


7.1 


1.4 


211.9 


NGC 7183 


SO 


-0.7 


13.14 


90 


close to box 




5 


77 


7200 


2.5 


0.5 


209.9 


NGC 7264 


Sb 


3.1 


14.58 


90 


elliptical 




2 


147 


13000 


1.2 


0.4 


248.0 


NGC 7332 


SO 


-2.0 


12.56 


90 


close to box 




1 


65 


15000 


7.5 


1.1 


121.3 


NGC 7396 


Sa 


1.0 


13.85 


60 


elliptical 




2 


13 


13000 


25.2 


1.6 


360.7 


NGC 7703 


SO 


-2.0 


14.41 


90 


elliptical 


cluster 


2 


57 


6400 


4.8 


2.0 


159.5 


NGC 7814 


SAab 


2.0 


11.68 


67 


elliptical 




1 


45 


18400 


33.5 


2.0 


231.0 


IC 1711 


Sb 


3.0 


14.53 


90 


elliptical 




2 


133 


6200 


2.1 


0.5 


173.4 


IC 1970 


SAb 


3.1 


13.36 


88 


elliptical 




4 


75 


7200 


8.1 


1.0 


118.9 


IC 2531 


Sb 


4.9 


12.88 


90 


peanut 


cluster 


3 


75 


7200 


2.3 


0.3 


228.1 


IC5176 


SABbc 


4.4 


13.02 


90 


elliptical 


companions? 


5 


151 


7200 


4.6 


0.3 


164.2 


IC 5264 


Sab 


2.3 


13.58 


90 


elliptical 


group 


4.5 


82 


12600 


2.8 


0.1 


230.6 


UGC 10043 


Sbc 


4.1 


15.17 


90 




cluster 


2 


61 


12200 


1.0 


0.3 


143.6 


UGC 11552 


Sab 


2.0 


15.68 


75 






5 


16 


7200 


11.7 


0.6 


243.0 


UGC 11587 


SO 


-1.9 


14.52 


90 






1 


76 


15000 


5.7 


0.7 




ESO 079-003 


SBb 


3.1 


13.91 


90 


close to box 




4.5 


134 


10800 


5.9 


2.0 


192.2 


ESO 234-053 


SO 


-1.9 


14.32 


90 






5 


81 


5400 


1.7 


0.1 




ESO 311-012 


SO-a 


0.1 


12.38 


90 


box 


isolated 


3 


14 


7200 


7.4 


1.3 


95.9 


ESO 443-042 


Sb 


3.1 


13.85 


90 


peanut 


companions? 


3 


130 


10800 


3.8 


0.7 


166.3 


ESO 512-012 


Sbc 


3.2 


14.42 


90 


peanut 




3 


113 


5400 


1.4 


0.3 


224.9 



a From Liitticke et al. (2000). 

b Run number in which each galaxy was observed. 

c Position angle of the long slit. 

d Total exposure time (in seconds). 

c Effective radius (in arcseconds) of the bulge, computed as explained in Section[4] 

f Maximum rotational velocity (in km s _1 ) along the disk major axis, corrected from inclination when necessary. 
Hubble type T, total apparent magnitude Bt, galaxy inclination i and Vmax have been taken from HyperLeda (http://leda.univ-lyonl.fr/ 1. 



The main parameters of the instrument configurations and 
other details of the observing runs are listed in Table [2] 

Using slit widths from 1.0 to 2.0 arcseconds, the different 
spectrographs allowed us to obtain spectra at a resolution of 
typically 4 A (FWHM, as measured from the widths of arc 
lines), except for run 3 (the ESO 3.6-m run), where a resolu- 
tion of 8.3 A was attained. 

The spectral range, from around A3900A to A5500A (with 
the exception of the ESO NTT runs, in which the spectra ex- 
tended redward to A6650A) was chosen to include the main 
prominent spectral features of the blue-yellow spectra of the 
galaxies, in particular most of the Lick/IDS line-strength in- 
dices. 

Table Q] includes some relevant observational parameters 
for the galaxy sample. The spectrograph slit was oriented along 
the minor axis of the bulges. 



Thus, effects of contamination by disk light are only ex- 
pected for the very central regions of the sample galaxies. 

Exposure times ranged from 1.5 to 5 hours per galaxy, pro- 
viding high signal-to-noise (S/N) spectra out to the outer re- 
gions of most bulges. 

We obtained reliable spectra out to the effective radius 
(r c ff) for every galaxy in our sample. Furthermore, our spec- 
tra reached S/N > 10 (per A and per arcsec in the spatial 
direction) beyond 2r c g for 80% of the sample galaxies. Some 
observations of nearby bulges were especially deep, providing 
useful spectra out to several kpc from the centres of the bulges 
(e.g., ~ 7 kpc for NGC 5084 and ~ 4 kpc for NGC 7814). 

In the central regions we obtained spectra with S/N ranging 
from 20 to 240, depending, among other factors, on the obscu- 
ration by dust lanes within the galaxy disk. 
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In each observing run we also acquired spectra of a num- 
ber of template stars from the Lick/IDS library (Worthey et al. 
1994) for calibration purposes (see below) as well as several 
spectrophotometric standards to calibrate the spectra in flux. 

4. Data reduction 

The reduction of the data was performed with the REDmlf] 
package (Cardiel 1999). We carried out a standard reduction 
procedure for spectroscopic data: bias and dark subtraction, 
cosmic ray cleaning, flat-fielding (using observations of tung- 
sten lamps and twilight skies to correct from high and low fre- 
quency variations respectively), C-distortion correction, wave- 
length calibration, S-distortion correction and recentering of 
the spectra, sky subtraction, atmospheric and interstellar ex- 
tinction correction, relative flux calibration, and spectrum ex- 
traction. RgDjnE is a package specifically written to reduce 
long-slit spectroscopic observation. Its main advantage is that, 
for each observed data frame (including calibration ones), it 
creates an associated error frame at the beginning of the re- 
duction procedure. From this point, error and data frames are 
processed in parallel, translating all the uncertainties in the ma- 
nipulation of the data into the error frames, and 

hence providing very accurate estimates of the random er- 
rors 

associated with the final spectra (see Cardiel et al. 1998). 

Although most of the reduction steps were performed us- 
ing standard procedures for long-slit reduction, we give below 
some comments on steps of particular importance: 

( i) Wavelength calibration. 

Spectra were converted to a linear wavelength scale using 
about 90 arc lines fitted by 5th to 7th-order polynomials. The 
rms dispersion of the fitted position of arc lines was typically 
0.25-0.30 A. 

( ii) Spectral rectification and centering. 

As usual with CCD spectrographs, the galaxy spectra were 
not completely perfectly aligned with the detector rows. To 
correct for this effect (which is crucial when measuring line- 
strength gradients), we fitted a Cauchy function to the number 
of counts in a spatial interval around the location of maximum 
counts for each spectral resolution element, hence determining 
the position of the galaxy centre as function of wavelength. In 
galaxies with dust lanes obscuring the light of the central re- 
gion of the bulge, these central parts were masked when fitting 
the luminosity profiles. Experiments with galaxies without dust 
lanes showed that we could still trace the position of the galaxy 
centre. The resulting map was then fitted by a low-order poly- 
nomial, which was in turn used to straighten the spectra. As 
a result of this process, the central spectrum of the galaxy is 
moved to the centre of a detector row so that symmetric spectra 
on each side of that row correspond to the same galactocentric 
radius. 

(Hi) Sky subtraction. 

Accurate sky subtraction is critical for studies like this 
since we aim to analyze spectra at light levels correspond- 
ing to only a few per cent of the sky signal. We refer the reader 



to Cardiel et al. (1995) for a description of the sky subtraction 
effects on the measurements of line-strength indices. For each 
of our galaxy observations we generated a sky image by fitting, 
for each pixel in the wavelength direction, a low-order poly- 
nomial using the regions close to the slit ends. In most cases, 
the galaxies fill only a small region of the slit, so this synthetic 
sky image is free from contamination from the galaxy light. 
However, this is not the case for some of the larger galaxies 
of the sample (like NGC 7814 and NGC 5084). In these cases, 
we fitted de Vaucouleurs' r 1 / 4 profiles to the surface brightness 
distribution of the galaxies along the slit, 

estimating the relative contribution of the galaxy to the 
regions from where the sky was extracted. We then 'decon- 
taminated' the sky frames by subtracting scaled and averaged 
galaxy spectra. In order to assess the quality of the sky subtrac- 
tion we looked for residuals due to sky emission lines, making 
sure that they were completely removed, even in the outer parts 
of the more extended galaxies. 

(iv) Measurement of effective radii. 

Effective radii for all the bulges of the sample were mea- 
sured by fitting de Vaucouleurs' profiles to the radial surface 
brightness profiles determined by collapsing the spectra in the 
wavelength direction. The central regions, affected in many 
cases by dust lanes, were removed from the fit. In any case, 
we never used spectra at radii below typically 2 arcsec. On the 
other hand, the fits were quite insensitive to the outer radius 
cutoff, although deviations of the luminosity profiles from a 
straight line in the logarithmic plots were sometimes observed. 
The derived effective radii are listed in Table Q] The associ- 
ated errors, in the last column, reflect the sensitivity of the 
derived radii to the particular choice of the region to be fit- 
ted. Note that we are not assuming that the bulges do follow 
a de Vaucouleurs luminosity profile (a Sersic profile with an 
exponent n closer to 1 provides a bett er fit in many cases, see 



Andredakis, Peletier, & Balcells 1995). However, the derived 



http://www.ucm.es/info/Astrof/software/reduceme/reduceme.html 



effective radii provide us with a well defined parameter to ex- 
press the line-strength gradients (Paper II) on an homogeneous 
spatial scale. 

(v) Flux calibration. 

Relative flux calibration of the spectra was achieved using 
exposures of Massey et al. (1988), Oke (1990) and Hamuy et 
al. (1994) spectrophotometric standards. In each run we took 
several repeated observations of 4, 5, 4, 3 and 3 stars (for runs 
1 to 5 respectively). In some cases we obtained exposures at 
different zenithal distances to check the atmospheric extinction 
correction. All the calibration curves of each run were averaged 
and the flux calibration errors were estimated by the differences 
among the indices measured with different curves. 

(vi) Extraction of spectra. 

Once the galaxy frames were fully reduced, we added all 
the frames of the same galaxies and extracted individual spec- 
tra along the slit. The availability of reliable error frames al- 
lowed us to compute the number of spatial bins to be coadded 
to guarantee a minimum S/N in the final spectra. Using the 
prescriptions given in Cardiel et al. (1998), we chose a lower 
S/N cutoff of 15 (per A; mean value along the full spectral 
range), which ensures a typical error in most of the atomic Lick 
indices below 20%. For the outer spatial regions, we coadded 
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Run 


Dates 


Telescope 


Instrument 


Dispersion 


Spectral range 


Slit width 


FWHM a 










A/pixel 


A 


arc sec 


A 


1 


13/08/96-16/08/96 


INT 2.5 m 


IDS 


1.6 


3908-5535 


2.0 


3.9 


2 


31/08/97-03/09/97 


INT 2.5 m 


IDS 


1.6 


3908-5535 


2.0 


4.3 


3 


08/03/97-11/03/97 


ESO 3.6 m 


EFOSC 


3.4 


3776-5504 


1.5 


8.3 


4 


23/09/97-25/09/97 


NTT 3.6 m 


EMMI 


1.3 


3971-6649 


1.0 


3.7 


5 


23/08/98-26/08/98 


NTT 3.6 m 


EMMI 


1.3 


3985-6649 


1.5 


4.9 



a Spectral resolution as measured from arc lines. 



some further spectra to get a minimum S/N of 10, which is 
enough to measure molecular indices (like Mg2) with a relative 
uncertainty below 10%, even though it does lead to significant 
errors in the atomic indices. 



5. Dynamical parameters 

Radial velocities and velocity dispersions for each spec- 
trum along the bulge radii were measured using a dedicated 
program within REDmE which incorporates the MOVEL 
and OPTEMA algorithms devised by Gonzalez (1993). The 
MOVEL procedure is a refined version of the classic Fourier 
quotient method (Sargent et al. 1977). The basic improvement 
over the standard procedure is that, starting from a first guess of 
7 (mean line-strength), V (radial velocity) and a (velocity dis- 
persion), it creates a model galaxy which is processed in par- 
allel to the galaxy spectrum. The differences between the input 
and output parameters of the model spectrum are then used to 
correct the galaxy power spectrum from systematic effects in 
the handling of data in the Fourier space. The OPTEMA algo- 
rithm allowed us to overcome the typical template mismatch 
problem. In order to do this, a number of stellar spectra of 
different spectral types and luminosity classes is fed into the 
program. In our case, we combined our template stars to con- 
struct representative input spectra of the following types: A3- 
5V, G0-2V, G5V, K0-1 V, G8-K0III, K3III, K5-7III, and M3III. 
The algorithm then constructs an optimal template as a linear 
combination of the above template spectra which minimizes 
the residuals between the galaxy spectrum and the broadened 
template. Meanwhile, the dynamical parameters are derived as 
part of the same minimization process. Our tests showed that if 
a single template star is used to measure the velocity dispersion 
(instead of the above combination), systematic errors as large 
as 10% can be introduced. 

Figure [2] illustrates the result of this procedure for two 
galaxies (NGC 7332 and UGC 10043), showing, in the up- 
per panels, the fits between the observed central spectra of 
the galaxies and the corresponding optimal templates corrected 
with the derived dynamical parameters. The lower panels show 
the residual of the fits. Note that emission lines, although al- 
most hidden in the original spectrum (like the [Olll]A5007 
line for NGC 7332) are clearly detected after subtracting the 
template (once an emission line is detected, the corresponding 
wavelength region is not used for the minimization of the resid- 
uals). In fact, the optimal templates derived for each galaxy 
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Fig. 2. Illustration of the velocity dispersion determinations 
using optimal templates for NGC 7332 (upper plot) and 
UGC 10043 (lower plot). For each plot, the upper panels show 
the observed galaxy spectrum (thick line) and the corrected op- 
timal template spectrum (thin line). The residuals of the fits are 
shown in the lower panels. 



were later used to mask the regions possibly containing emis- 
sion lines. 

Realistic errors in the derived parameters (radial velocity 
and velocity dispersion) were computed by performing Monte- 
Carlo simulations, repeating the whole process (including the 
derivation of the optimal template) for a large number of simu- 
lated galaxy spectra created using the error spectrum obtained 
during the reduction process. We checked that, in order to mea- 
sure velocity dispersion of the order of the spectral resolution 
(~ 100 km s _1 for runs 1, 2, 4 and 5) with a relative error be- 
low 20% we needed spectra with S/N (per A) larger than 15. 
This coincides with the cutoff value for the spectral extraction 
extraction to measure the atomic indices. 
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Once an optimal template is derived for the central or max- 
imum S/N (if the centre is affected by dust lanes) spectrum, 
the procedure is repeated for all the spectra along the bulge 
radii, deriving radial profiles for the radial velocity and veloc- 
ity dispersion. Note that the template derived for the central 
spectrum is used for the analysis of the rest of the spectra and 
that, therefore, that we are not taking into account the possi- 
ble effect of population gradients in the derived sigma profiles. 
However, we preferred to introduce this small inaccuracy than 
to artificially increase the scatter of the parameters profiles due 
to unreal changes in the derived linear combinations of tem- 
plates at each radius. Note also that any radial change in mean 
line-strength (i.e. 7) or continuum shape is properly corrected 
from in the MOVEL algorithm. 

The dynamical parameters, and their corresponding errors, 
at each galactocentric radius are listed in Table [3] Since the 
slit was oriented along the minor axes of the bulges, no ro- 
tation along the slit was expected. To check this, in Table [3] 
we include the maximum rotational velocity (computed as an 
error-weighted mean of the, generally two, spectra at each side 
of the galaxy centre with the highest rotational velocity) and its 
error. It is apparent from the Table that most bulges do not show 
any significant rotation along the minor axes. For some galax- 
ies, however, we observed some hints of rotation or peculiar V r 
profiles in the central regions, probably due to slight misalign- 
ments of the slit with respect to the minor axis direction or the 
bulge centre. 

For galaxies with no obvious dust lanes in the central re- 
gions, the table also includes the central velocity dispersion <t > 
measured in an aperture of 4 arcsec times the slit width centred 
on the galaxy centre. We also list mean velocity dispersions 
((a)) defined as the error- weighted mean velocity dispersion of 
all the extracted spectra along the radius (which is very close 
to the velocity dispersion that is derived when adding all the 
spectra along the slit). 

The derived velocity dispersion profiles for the galaxy sam- 
ple are plotted in Figure[3] It can be seen that, typically, the ve- 
locity dispersion decreases with galactocentric radius, although 
there are different behaviors and some galaxies exhibit a large 
scatter. It should be noted that the main motivation for mea- 
suring velocity dispersion profiles in this paper is to have the 
necessary information to carry out accurate corrections to the 
Lick/IDS indices (detailed further below). 

Since line-strength indices depend on the spectral resolu- 
tion, all spectra along the slit must be corrected to the same 
velocity dispersion before measuring the indices so as to avoid 
the introduction of spurious line-strength gradients. However, 
if we used the nominal velocity dispersion values obtained for 
each individual spectrum to perform this correction, the errors 
associated with the a measurement would be introduced in the 
measured indices. 

To avoid this, we smoothed the velocity dispersion pro- 
files by fitting low-order polynomials to the data, and we used 
the predictions of these polynomials for correcting the indices 
to the desired velocity dispersions. Note that we are not pre- 
tending to assign any physical meaning to these polynomial 
radial profiles; they only have a calibration purpose. The col- 



umn er max in Table [3] gives the maximum velocity dispersion 
indicated by the polynomials in the observed radial range. 

Finally, we compared the derived values with the avail- 
able data from the literature (last columns of Table 0) to check 
for any systematic offset in the derived velocity dispersion. To 
minimize aperture differences, these reference data should be 
compared with our central (ctq) values when available, or with 
(a) otherwise. It can be seen that we do not have any systematic 
offset relative to the scarce previous data. 



6. Line-strength indices 

Line-strength indices in the Lick/IDS system (e.g. Worthey et 
al. 1994; Worthey & Ottaviani 1997) were measured for all the 
extracted spectra along the measured extent of the bulge. We 
measured all the indices from H8 to Fe5406. In the case of 
runs 4 and 5 (cf. Table |2), the extended spectral range allowed 
us to include also the redder Lick indices (up to the TiO in- 
dices). The next subsections explain all the procedures to derive 
reliable errors in the line-strengths and to transform them to 
the spectrophotometric Lick/IDS system. All these steps were 
possible thanks to the inclusion of a number of stars from the 
Lick/IDS library in the different observing runs. In particular, 
we observed 40, 40, 19, 4 and 6 Lick standard stars (covering 
a range of spectral types and luminosity classes) in runs from 1 
to 5 respectively. 

6. 1. Error analysis 

A first estimation of errors in the line strengths was accom- 
plished by taking into account three different error sources: 
(i) photon noise (the errors were calculated using the vari- 
ance spectra computed as part of the reduction process), (ii) 
flux calibration (the estimation was based on the comparison of 
the indices measured using the different flux calibration curves 
obtained in each observing run), and (//^uncertainties in the 
wavelength scale (due to errors in the wavelength calibration 
and in the measurement of the radial velocity). 

Except for the spectra with the highest S/N, the dominant 
error source was the first one. 

In order to check and refine the error analysis, we ana- 
lyzed the line-strengths obtained for: ( i) repeated observations 
of the same stars within each observing run, ( ii) stars in com- 
mon among different observing runs (e.g. all the 40 Lick stars 
observed in run 1 were repeated in run 2), and (Hi) separate 
observations of the same galaxies within each run (typically 
the observation of a galaxy was splitted into 4-12 different ex- 
posures). Using this data, we computed residual errors as the 
quadratic differences between the measured r.m.s. dispersions 
and the expected total error from the three sources of uncer- 
tainty quoted in the previous paragraph. The resulting residual 
errors are listed in Table [5] for each index and observing run. 
They are relatively modest but were included anyway in the 
final uncertainties. 
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Fig. 3. Velocity dispersion profiles for the galaxy sample. Filled and open symbols refer to spectra of the bulges at either side of 
the galaxy centres. The plotted lines represent polynomial fits to the profiles. See text for details. 
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Table 3. Dynamical parameters of the bulge sample. 



Galaxy 


K a 


AVi 


r 

V 4- 

v rot 


AVrot 




Acr 


Ti 

< o > d 


A < a > 


e 

'-'max 


G ref 


T— 

Ref 


NGC 522 


2719.0 


12.5 


12.1 


4.9 


82.1 


9.3 


88.2 


8.8 


99 






NGC 585 


5358.4 


6.8 


5.5 


4.5 






174.3 


4.7 


188 






NGC 678 


2814.9 


15.0 


13.3 


6.3 


163.4 


2.1 


170.7 


4.1 


183 






NGC 891 


548.2 


12.3 


1.9 


2.2 






95.4 


5.1 


95 


72 


(2) 


NGC 973 


4081.0 


9.6 


1.5 


6.0 






173.3 


9.3 


226 






NGC 1032 


2702.9 


7.7 


16.3 


9.4 






214.2 


2.8 


225 






NGC 1184 


2278.2 


8.0 


5.1 


3.1 


230.3 


3.0 


219.6 


7.6 


234 


229 


(3) 


NGC 1351A 


1329.3 


17.6 


6.1 


15.4 






141.4 


28.5 


141 






NGC 1886 


1689.8 


15.1 


8.0 


3.3 






111.0 


6.5 


120 






NGC 3957 


1589.7 


10.2 


1.6 


2.2 






131.4 


5.3 


149 






NGC 5084 


1671.6 


16.8 


10.5 


9.5 


207.6 


1.7 


200.0 


4.0 


214 


215 


(1) 


NGC 6010 


2058.5 


3.3 


2.9 


3.3 


148.4 


1.2 


146.7 


4.9 


157 


144 


(4) 


NGC 6829 


3303.3 


11.0 


1.2 


5.6 






118.1 


6.3 


130 






NGC 7183 


2594.9 


12.8 


2.2 


5.0 






138.8 


5.3 


157 






NGC 7264 


4210.1 


14.4 


6.9 


4.6 






165.6 


6.9 


205 






NGC 7332 


1206.3 


2.2 


2.6 


2.8 


117.2 


0.6 


116.1 


2.8 


122 


134 


(1) 


NGC 7396 


4882.5 


12.2 


8.9 


6.6 






220.6 


4.1 


224 






NGC 7703 


3970.5 


7.1 


7.1 


4.8 


149.9 


3.2 


148.5 


5.9 


154 


166 


(3) 


NGC 7814 


1051.4 


6.9 


4.8 


4.5 






170.3 


2.4 


187 


172 


(1) 


IC 1711 


2785.2 


11.8 


9.0 


5.2 






122.3 


12.2 


161 






IC 1970 


1300.8 


7.2 


11.6 


10.6 






82.8 


12.1 


95 






IC 2531 


2386.8 


11.7 


3.8 


4.1 






126.6 


8.7 


127 






IC5176 


1739.3 


10.5 


4.1 


4.3 






152.7 


6.4 


160 






IC 5264 


2025.1 


6.3 


2.6 


3.6 






101.1 


3.4 


106 






UGC 10043 


2180.3 


16.0 


8.2 


6.4 






88.5 


15.7 


89 






UGC 11552 


4479.5 


11.1 


2.7 


11.3 






151.3 


10.7 


167 






UGC 11587 


4481.5 


7.8 


2.3 


4.9 


213.0 


2.9 


203.4 


8.5 


224 






ESO 079-003 


2613.5 


8.1 


1.5 


6.8 






120.7 


8.5 


160 






ESO 234-053 


6232.9 


16.0 


0.3 


4.9 


204.9 


3.1 


200.4 


7.6 


226 






ESO 311-012 


1106.8 


10.7 


3.2 


5.3 


124.2 


3.5 


132.1 


4.0 


147 






ESO 443-042 


2815.8 


24.8 


8.7 


9.6 






118.6 


14.9 


119 






ESO 512-012 


3313.4 


19.0 


13.0 


6.5 






122.8 


12.4 


123 







a Mean radial velocity (in km s _1 ). 

b Maximum rotational velocity (see the text). 

c Central velocity dispersion. 

d Mean velocity dispersion. 

c Maximum velocity dispersion (see the text). 

f References are: (1) Average value quoted in HYPERLEDA; (2) Bottema et al. (1991); (3) Di Nella et al. (1995), and (4) Falcon-Barroso 
et al. (2002). 



6.2. Emission-line corrections 

Some line-strength indices, in particular H/3, Fe5015 and Mgb, 
are potentially affected by emission lines (e.g., Goudfrooij & 
Emsellem 1996). Most of our galaxies exhibit non-negligible 
nebular emission, as judged from their [O III] A5007 line. 

Usually, the emission is concentrated towards the bulge 
centres, where the spectra are heavily contaminated by the 
disk light (Hll regions can project onto the slit). However, 
in some special cases (e.g., NGC 1886 and NGC 7332), 
line emission can be detected even in the outer part of the 
bulge, where any contamination by disk light is negligible. 
For instance, in NGC 7332, the line emission increases to 
£W([Oill]A5007) ~ 4 A at a distance of - 16 arcsec from 
the centre. 



The correction of H/3 from the emission contribution is a 
troublesome task. The usual procedure (see eg. Gonzalez 1993, 
Trager et al. 2000) relies on using the equivalent width of [O III] 
to apply an empirical correction. In this paper, given the un- 
certainties of that correction and the fact that the emission is 
relatively strong for many of the galaxies, we have decided not 
to correct H/3 and, therefore, not to use it in our analysis in the 
cases where any [Om]A5007 emission was detected (in par- 
ticular when its measured equivalent width was larger than its 
error). 

We determined equivalent widths of [O III] A5007 after sub- 
tracting from each galaxy spectrum an emission-free template 
spectrum constructed from the broadened optimal template de- 
rived for that galaxy during the velocity dispersion measure- 
ment (see Figure|2|i. 
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In order to measure the Fe5015 and Mgb indices in spectra 
with emission lines, we masked the spectral regions affected 
by the [Om]A4959, [Olll]A5007, and [Nl]A5 199 lines and re- 
placed them with the corresponding regions in the broadened 
optimal template spectrum. 

6.3. Broadening corrections 

As mentioned before in Section|5] the Lick/IDS indices have to 
be corrected for spectral broadening. In order to avoid 

the introduction of artificial line-strength gradients due to 
changes in velocity dispersion along the slit, we broadened all 
the spectra of a given galaxy to match the maximum velocity 
dispersion <r max derived in Section [5] by convolving it with a 
Gaussian of width (o^ lax — of) 1 / 2 , where ay is the velocity 
dispersion at each radius predicted by the corresponding poly- 
nomial fit. 

The conversion of the line indices to the resolution of the 
Lick system (see Section Rvfl i was done as follows. In the cases 
where the total resolution of the spectra (the quadratic sum 
of instrumental resolution and the maximum velocity disper- 
sion for the galaxy) was below the Lick resolution, we sim- 
ply broadened the spectra by the quadratic difference. In other 
case, we applied an empirical correction to the measured in- 
dices. This correction was derived by broadening stellar spec- 
tra by convolution with Gaussians having a range of a values. 
Since this correction depends on the individual spectral type of 
the star, we used instead the optimal templates obtained in the 
velocity dispersion measurement procedure to derive the fol- 
lowing expression: 

I (era) = 1(a) (ci + c 2 a/a Q + c 3 (a/a ) 2 ^ , (1) 

which can be used to convert an index / measured at a resolu- 
tion a to any other resolution erg. The values of the Cj coeffi- 
cients are listed in Table|4]for all the Lick indices. In principle, 
the exact values of the coefficients depend on the particular op- 
timal template used. However, the spectral differences among 
them have a small effect on the derived polynomials and the 
uncertainties introduced by using the listed mean coefficients 
are well below the other sources of errors. 

6.4. Conversion to the Lick/IDS system 

Transformation of line-strengths to the Lick/IDS system is not 
straightforward. Basically, there are two important effects that 
have to be taken into account. First, the Lick/IDS spectra are 
not flux-calibrated using a spectrophotometric standard star, 
but by using a normalized tungsten lamp spectrum. Thus, there 
are significant differences in the spectral shapes when com- 
pared with flux-calibrated spectra. This effect causes system- 
atic offsets, especially for indices with broad bands like (e.g.) 
Mg2. Second, the spectral resolution of the Lick spectra is not 
constant with wavelength, in particular it degrades consider- 
ably towards the blue end. This has an important effect on 
narrow line indices, like many of the atomic ones. For a com- 
prehensive discussion on these effects, we refer the reader to 
Worthey & Ottaviani (1997). 



Table 4. Coefficients of the polynomials to correct Lick indices 
for broadening effects. 



Index 


Cl 


C2 


C3 




0.76129 


0.10120 


0.13751 


H5 f 


1.25993 


-0.51501 


0.25508 


CNi 


0.70442 


0.31941 


-0.02382 


CN 2 


0.76492 


-0.00441 


0.23949 


Ca4227 


0.78022 


-0.28090 


0.50067 


G4300 


0.89152 


0.10002 


0.00847 


H 7A 


1.11845 


-0.21248 


0.09403 


H7F 


0.80429 


0.25651 


-0.06080 


Fe4383 


0.87985 


-0.00952 


0.12967 


Ca4455 


0.63849 


0.09619 


0.26532 


Fe4531 


0.91332 


-0.01235 


0.09903 


Fe4668 


1.00545 


-0.12685 


0.12140 


H/3 


0.99068 


0.00563 


0.00369 


Fe5015 


0.78157 


0.17915 


0.03928 


Mgi 


0.97728 


0.01066 


0.01206 


Mg 2 


0.98082 


0.02199 


-0.00282 


Mgb 


0.96450 


-0.07495 


0.11045 


Fe5270 


0.82533 


0.12283 


0.05184 


Fe5335 


0.84324 


-0.08141 


0.23817 


Fe5406 


0.89017 


-0.14121 


0.25104 


Fe5709 


0.89437 


0.01365 


0.09198 


Fe5782 


0.86475 


-0.14458 


0.27983 


Na5895 


0.94445 


0.01617 


0.03938 


TiOi 


1.00000 


0.00000 


0.00000 


T1O2 


1.00000 


0.00000 


0.00000 



To transform our data to the Lick/IDS system we performed 
the following steps. First, we estimated the resolution at which 
each particular index should be measured by broadening our 
large set of stars in common with the Lick library to several 
line widths, in steps of 25 km s _1 . We did not attempt to find 
precise resolutions for each index since these depend on the 
particular set of stars being used. Instead, we estimated the ap- 
proximate line widths ol which, changing softly with wave- 
length, minimized the residuals. These are given, for each in- 
dex, in the second column of Table [5] These resolutions agree 
with the rough estimates given in Worthey & Ottaviani (1997). 

We then broadened all our stellar spectra to the above 
<jl resolutions and measured the line-strength indices. The 
comparison of these measurements with those in the original 
Lick/IDS spectra allowed us to derive mean offsets for all the 
indices in each observing run. These additive offsets, listed in 
Table [5] are mainly due to differences in flux calibration be- 
tween both systems, but also include a fine tuning of the broad- 
ening corrections and systematic effects in both our data (like 
uncertainties in the spectral resolution or flux calibration) and 
the Lick/IDS data (like spectral peculiarities which depend on 
the particular Lick/IDS run in which each star was observed). 
Note that we obtain different offsets for each observing run (al- 
though in most cases they are similar). In fact, the comparison 
with Lick data provides us with an anchor system to correct for 
small offsets among the different runs. It must be noted that the 
offsets for run 3, with a significantly lower spectral resolution, 
are of the same order than the offsets for the other runs, which 
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confirms that the broadening corrections have been properly 
done. 

Figure|4]shows a comparison of the line-strength indices in 
the Lick/IDS spectra and the measurements in our sample of 
stars after correcting to the Lick/IDS system. Error bars where 
derived by combining the computed errors for our measure- 
ments and the errors of the original Lick/IDS indices as taken 
from Worthey et al. (1994). The last two columns of Table [5] 
give the r.m.s. standard deviation with respect to the zero off- 
set, and the standard deviation that should be expected from the 
computed errors in the offsets. As it can be seen, the former are 
usually larger than the latter, indicating that the errors in the 
/(Lick) — /(here) differences are somewhat underestimated. 
Given all the tests that we have carried out with our data, we 
think that this disagreement is most likely due to an underesti- 
mate of the errors in the original Lick measurements. 

6.5. Final data 

Once all the previous corrections were derived, we applied 
them to the extracted spectra along the bulge radii and mea- 
sured the line strengths. We derived two sets of final data: 

1 . In the first one, we introduced the residual errors, masked 
the emission lines, converted all the indices to a resolution 
corresponding to the maximum velocity dispersion along 
the radius (as listed in Table Up, and applied relative off- 
sets among runs to put them in the same spectrophotometric 
system (which should be very close to a true flux-calibrated 
system). To do this, we applied the offsets listed in Table 
5 with respect to run 2 (offsets for run 2 were therefore 
zero). However, we did not transform the indices to the 
Lick/IDS system. This data can be very useful for future 
studies working at any spectral resolution and in a flux- 
calibrated system. Note that, in this dataset, indices from 
different galaxies cannot be directly compared, since in- 
dices are measured at different resolutions. 

2. For the second dataset, we carried out a full transformation 
to the Lick/IDS system, applying (in addition to the correc- 
tions from the last paragraph) the necessary offsets and the 
index-dependent velocity dispersion corrections. For this 
last step, we broadened the spectra when possible or ap- 
plied the polynomials of Table [4] when the original resolu- 
tions (the quadratic sum of cr max and the instrumental res- 
olution) were larger than ox (Table[5]). 

6.6. Comparison with other works 

Finally, in Figure [5] we present a comparison of the line- 
strengths measured in the central regions of our bulges with 
the available results from the literature. The galaxies with ex- 
isting published data and the corresponding references are the 
following: NGC 5084 (Trager et al. 1998), NGC 6010 (Falcon- 
Barroso et al. 2002), NGC 7703 (Caldwell, Rose & Concannon 
2003), NGC 7332 (Bender, Burstein & Faber 1993; Trager et 
al. 1998; Golev et al. 1999; Falcon-Barroso et al. 2002, 2004; 
Denicolo et al. 2005; Sil'chenko 2006; Kuntschner et al. 2006; 
Sanchez-Blazquez et al. 2006), and NGC 7814 (Prugniel et al. 



2001). To perform the comparison, for each galaxy and refer- 
ence we binned our spectra in the central regions of the bulges 
to reproduce as closely as possible the aperture used by the 
corresponding works. It is apparent from the figure that, al- 
though there is a general agreement with the data from other 
workers, there are also some important discrepancies for some 
particular line-strengths and authors (like some indices for 
NGC 5084 from Trager et al. 1998, or the Mg2 measurement 
for NGC 7814). These differences could be due to: low S/N ra- 
tios in the spectra of other authors, variations in the centering of 
the slit or in the observed aperture, or differences in the calibra- 
tion and data reduction processes. However, the fact that, when 
there are several published values for a galaxy (NGC 7332), our 
measurements agree with the majority of them gives us confi- 
dence in our measured values. For instance, for the 5 previ- 
ous measurements for the central Mg2 index in NGC 7332, our 
value agrees within the errors with 4 out of them. 

The comparison with the central line-strengths from 
Sanchez-Blazquez et al. (2006) is of particular relevance, since 
in Paper II we will make use of this dataset as a comparison 
sample of elliptical galaxies. In this case, in which we plot 
with blue symbols the indices measured in an aperture of 2 x 4 
arcsec (as listed in Sanchez-Blazquez 2004), we find a good 
agreement with our data (with the exception of the H7F index, 
although note that the results agree for H7A). 



7. Summary 

We have introduced our long-term project devoted to under- 
stand the evolutionary status of galaxy bulges by presenting the 
sample and the most relevant details of our survey of long-slit 
spectra along the minor axis of bulges in edge-on spiral galax- 
ies. After summarizing a data reduction process characterized 
by a very detailed error analysis, we have presented our pro- 
cedures to measure dynamical parameters and all the Lick/IDS 
line-strength indices available in the spectra. A special empha- 
sis has been made in the broadening corrections of the indices 
and their conversion to the Lick/IDS system. As a final product, 
we have derived tables giving all the parameters (radial veloci- 
ties, velocity dispersions and indices) at each galactocentric ra- 
dius for all the bulges of the sample. In the case of the indices, 
two datasets have been obtained: one with the indices in flux- 
calibrated spectra and at the natural resolution of the bulges 
(i.e. internal plus instrumental), and another in which all the 
indices have been converted to the Lick/IDS system. This two 
datasets will be available in electronic form at CDS. 

In the next paper of this series we will measure line- 
strength gradients and analyze them at the light of different 
galaxy formation scenarios. 
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Table 5. Spectral resolutions, residual errors and offsets needed to transform line-strength indices to the Lick/IDS system. 



Index 


°l 
kms" 1 


run 1 


Residual errors 
run 2 run 3 


runs 4,5 


run 1 


run 2 


Offsets 
run 3 


runs 4,5 


rms 


exp 


H<5 A 


325 


0.00 


0.00 


0.00 


0.20 


-0.34 


-0.34 


0.00 


-1.09 


0.66 


0.37 


FMf 


325 


0.00 


0.00 


0.00 


0.00 


-0.12 


-0.09 


0.00 


-0.81 


0.35 


0.24 


CNi 


325 


0.005 


0.005 


0.005 


0.008 


+0.003 


+0.003 


-0.009 


0.000 


0.020 


0.015 


CN 2 


325 


0.004 


0.004 


0.004 


0.011 


+0.008 


+0.008 


-0.015 


-0.011 


0.021 


0.016 


Ca4227 


300 


0.04 


0.04 


0.04 


0.00 


0.00 


0.00 


-0.08 


0.00 


0.23 


0.17 


G4300 


300 


0.13 


0.13 


0.13 


0.00 


-0.08 


-0.08 


0.00 


0.00 


0.38 


0.32 


H 7A 


275 


0.19 


0.19 


0.21 


0.00 


+0.41 


+0.41 


+0.41 


+0.34 


0.48 


0.40 


H7F 


275 


0.04 


0.00 


0.00 


0.00 


+0.17 


+0.17 


0.00 


+0.10 


0.22 


0.19 


Fe4383 


250 


0.22 


0.21 


0.21 


0.09 


+0.08 


+0.08 


0.00 


-0.36 


0.53 


0.45 


Ca4455 


250 


0.04 


0.04 


0.04 


0.00 


+0.32 


+0.32 


0.00 


+0.28 


0.20 


0.20 


Fe4531 


250 


0.10 


0.10 


0.10 


0.18 


+0.15 


+0.15 


-0.32 


0.00 


0.35 


0.31 


Fe4668 


250 


0.15 


0.19 


0.15 


0.57 


0.00 


0.00 


-1.03 


0.00 


0.57 


0.59 


H/3 


225 


0.10 


0.10 


0.10 


0.04 


-0.07 


-0.07 


0.00 


-0.14 


0.22 


0.20 


Fe5015 


200 


0.15 


0.19 


0.15 


0.21 


+0.26 


+0.26 


-0.11 


0.00 


0.64 


0.40 


Mgi 


200 


0.004 


0.004 


0.004 


0.007 


+0.011 


+0.003 


+0.016 


+0.009 


0.011 


0.010 


Mg 2 


200 


0.006 


0.005 


0.003 


0.005 


+0.025 


+0.020 


+0.021 


+0.016 


0.012 


0.010 


Mgb 


200 


0.04 


0.03 


0.02 


0.07 


0.00 


0.00 


-0.08 


0.00 


0.31 


0.17 


Fe5270 


200 


0.00 


0.00 


0.00 


0.11 


0.00 


+0.03 


0.00 


0.00 


0.23 


0.16 


Fe5335 


200 


0.02 


0.02 


0.02 


0.16 


-0.06 


-0.06 


0.00 


+0.18 


0.26 


0.16 


Fe5406 


200 


0.04 


0.04 


0.04 


0.03 


-0.04 


-0.04 


-0.19 


0.00 


0.23 


0.14 


Fe5709 


200 








0.00 








-0.08 


0.14 


0.14 


Fe5782 


200 








0.00 








+0.14 


0.10 


0.13 


Na5895 


200 








0.24 








0.00 


0.26 


0.40 


TiOi 


200 








0.003 








+0.013 


0.008 


0.007 


TiQ 2 


200 








0.000 








-0.015 


0.005 


0.004 



a Spectral resolutions at which line-strength indices should be measured to roughly match the original Lick/IDS system. 

b Measured (rms) and expected (exp) standard deviations for the comparison of indices with the original values of the Lick/IDS system. 
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Fig. 4. Comparison of the line-strength indices in the Lick/IDS spectra and the measurements in our sample of stars after correct- 
ing to the Lick/IDS system. Different symbols refer to the different observing runs, using the following code: circles for runs 1 
and 2; triangles for run 3, and stars for runs 4 and 5. In the case of the TiO indices, both TiOi are T1O2 are plotted in the same 
panel, using closed and open symbols for both indices respectively. The length of the error bar near the left end of each panel is 
twice the r.m.s. standard deviation of each dataset with respect to the zero offset horizontal line. 




Fig. 4. (Continued) 
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